6 Li in the atmosphere of GJ 117 
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^ ■ ABSTRACT 



We present high resolution VLT UVES observations of the active K dwarf 
GJ 117. 6 Li enhancement has been shown for energetic solar events, one chromo- 
spherically active binary, and several dwarf halo stars. Our analysis reveals the 



detection of 6 Li on this source with = 0.030±0.010. We found no significant 
contribution from other lines, including Ti I, in the Li profile of GJ 117 and a 
template star of similar spectral type and metallicity. We discuss the possibility 
for 6 Li production by spallation and find it to be consistent with the activity 
levels of the object. 



Subject headings: Stars: activity - Stars: atmospheres - Stars: individual: GJ 
117 -Stars: late-type 
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1. Introduction 

Lithium is a fragile element that serves as a powerful diagnostic of stellar structure and 
evolution. The fragility of Li arises from the fact that it can burn at relatively low tem- 
peratures via the [ 7 Li (p,a) 4 He] and [ 6 Li (p,a) 3 He] processes. Most of the Li destruction 
occurs during the pre-main sequence state, when the convection zones are deepest and for 
stellar models with the lowest mass. Lower mass stars arrive on the ZAMS at lower temper- 
atures and the amount of protostellar Li that survives today in dwarf stars decreases with 
decreasing temperature (Deliyannis & Malaney 1995). Stellar clusters such as a Per and the 
Pleiades, show a large spread of Li abundances at each effective temperature and it has been 
shown that the largest Li abundances are found in fast rotating stars with the strongest chro- 
mospheric emission (Soderblom et al. 1993). This is a puzzling result as theoretical models 
predict that fast rotation enhances the mixing process that leads to increased Li depletion; 
the rotational history plays a major role in the depletion process (Pinsonneault 1997). 

On the main sequence, the correlation between Li and chromospheric activity is not 
very strong. Stars with the same Ca II flux may have Li abundances that differ by one 
order of magnitude. Lithium is enhanced in sunspots (Giampapa 1984) where a small but 
notable amount of 6 Li has been detected (Ritzenhoff et al. 1997). An investigation of active 
RS CVn binaries failed to show a clear correlation between photometric variations and Li 
abundances (Pallavicini et al. 1992). Plages and active regions could play a minor role 
but their effect would be to decrease rather than increase the Li 6707.8A equivalent width 
(Houdebine & Doyle 1995). Russell (1996) used the Li 6104A line, which is formed deeper 
into the atmosphere, to study a sample of Pleiades K dwarfs. The results showed that the 
derived abundances have significantly less scatter than those derived from the 6707. 8A line. 
If activity has an effect on Li, it seems to increase rather than decrease the equivalent width 
and is therefore opposite to what the chromospheric NLTE calculations predict. 

Surface activity could increase the Li abundance if extreme energetic conditions are 
met (see Mullan & Linsky 1999). Spallation reactions from low energy flares can increase 
ratios (Canal et al. 1975). Predictions for significant creation of 6 Li in stellar and solar 
flares have been made (Canal et al. 1975; Walker et al. 1985; Livshits 1997) although actual 
measurements for the Sun have found the average rj^ °f 0-01 and up to 0.04 (Traub & Roesler 
1971; Miiller et al. 1975). 

The first detection of 6 Li on the halo dwarf HD 84937 has been highly significant (Smith 
et al. 1993). Follow up studies have confirmed this detection and have also shown the presence 
of the isotope in BD+26 3578 (Smith et al. 1998; Hobbs & Thorburn 1997). Unexpected 
6 Li abundances have been observed during a long flare on a chromospherically active binary 
(Montes & Ramsey 1998). These have been attributed to spallation reactions. 
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In this paper we present high resolution VLT UVES observations of the K star GJ 117 
with emphasis on the Li 6707.8A line profile. GJ 117 (=HD 17925) is a single chromospher- 
ically active K dwarf. The object shows evidence for photometric variability with a period 
of 6.5 days (Cutispoto 1992) while the weak color variations show that the star is redder 
at minimum light. Spectroscopic observations in the Ca II lines also show strong evidence 
for variability at both short and long timescales (Baliunas et al. 1995). Our observations 
and analysis are presented in § 2. In § 2.1 we present stellar models generated using the 
PHOENIX code for the temperature and gravity of GJ 117. We present the results, derived 
lithium parameters, and ^ ratio in § 2.2. Discussion of the results and implications for 
lithium production are presented in § 2.3, and lastly we summarize our findings in § 3. 



2. Observations & Data Analysis 

Measurements of 6 Li in stellar spectra are very difficult and rely on the modeling of 
the increased width and red asymmetry of the 6707.8A doublet. The problem becomes 
particularly difficult in active late-type stars where large rotational velocities will smear the 
line profiles. The problem is further complicated by the occurrence of blends, such as Ti I 
and CN. The observations were conducted in December 2002 with the VLT Kueyen Telescope 
and UV- Visual Echelle Spectrograph (UVES). This set-up used the 600 1 mm -1 red grating 
and MIT/LL eev 2kx4k CCD. With an exposure time of 150 seconds this setup provided a 
resolution of ~1 10,000 and a signal-to- noise of f«400 in the final spectrum. The VLT data 
were reduced with the VLT data pipeline (uves/1.3.3) with further analysis carried out using 
the STARLINK based DIPSO software. The UVES spectrum of GJ 117 in the Li I 6708A 
region is shown in Figure 1. 



2.1. Model Spectra 

Theoretical spectra were calculated using the general stellar atmosphere code Phoenix 
(see Hauschildt et al. 1995; Allard & Hauschildt 1995). We have used the LTE atmosphere 
models from the NextGen series of (Hauschildt, Allard, & Baron 1999) for the effective 
temperature and gravity listed in Table 1. Given the uncertainties in T e ff and \og(g), 
typically ± 200K and ±0.3 dex, we have used the directly converged models from NextGen 
which are sampled every 200K in T e ff and 0.5dex in \og(g). Line blanketing is handled 
within Phoenix by means of direct opacity sampling (Hauschildt et al. 2001). 

The 6 Li resonance doublet has been included in Phoenix by adding the wavelengths and 
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gf values of Smith et al. (1993) into the master line list Kurucz (1995) where the fraction of 
is introduced as a parameter (7). In direct opacity sampling, opacity at each wavelength 
point is calculated as a sum of opacities from all the contributing species. The contributions 
from 7 Li and 6 Li are then expressed as : 

k( 7 Li) = {\- 1 ) K{ 7 Li tot ) (1) 

k( 6 Li) = 7 K( 6 Li tot ) (2) 

where k( x Li tot ) is the opacity calculated as if all lithium is in the X isotopic state. Theoretical 
line profiles have been calculated with 7 ranging from 0.0 to 0.2. Microturbulence has been 
taken into account when generating the model spectra using a velocity of £=1.5 km s _1 . The 
v sini of the object was derived from x 2 analysis of several lines in the Li I 6707. 8A region. 

2.2. The ^ isotope ratio 

We have carried out a comparison between the observed and model profiles using % 2 
statistics. A grid of models in steps of yj^ = 0.01 was developed and the most probable 
model with x 2 ~1 was determined. Figure 2 shows the Ax 2 computed for the model grid 
and plotted with the sign preserved. The best fit has a £g = 0.030 ± 0.007 and Ax 2 ~ 30 
as compared to a model with no 6 Li. The observed spectrum of GJ 117 fitted with ^4 ratios 
of 0.0, 0.03 and 0.10 is shown in Figure 3. x 2 was calculated for 35 degrees of freedom and 
a sign factor was added to x 2 to indicate whether the observational data was larger than 
(positive) or smaller than (negative) the model. 

The ratio derived from the best fits are also given in Table 1. la errors for the 
table were derived using x 2 + 1 e.g. Nissen et al. (1999), and we also include a conservative 
5% systematic error to account for the uncertainty in the continuum level making 7^ = 
0.030±0.010. 

The influence of lines such as Fe I, Ti I, and CN on the Li I profile, can potentially 
mimic an enhanced 6 Li abundance (Nissen et al. 1999; Reddy et al. 2002). The effects of 
possible CN bands near 6707.5 have been seen to be < 1 mA in the solar spectrum and 
should have a negligible effect for GJ 117. In earlier lower resolution data the Fe I line at 
6707.4 A was often blended with Li I, however in the UVES observations these are easily 
separable. We have tested models with enhanced Fe abundances and found that these have 
no effect on the Li I line profile. 
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The most significant blend may come from the Ti I lines in the vicinity of 6708A and 
could lead to erroneous results. Israelian et al. (2001) have claimed the detection of ^ = 
0.12 in the planet hosting star HD82943. Reddy et al. (2002) concluded that this enhanced 
6 Li could be explained by the presence of a Ti I lines near 6708.03 and 6708.1 A. In a 
follow up paper Israelian et al. (2003) revised the ratio to £g = 0.05±0.02. The PHOENIX 
models used in our analysis were updated to include the atomic data for Ti I lines taken 
from Table 3 of Reddy et al. (2002). We have investigated the effects of Ti I on the Li line 
profile by constructing a grid of models with Ti abundances ranging from log(Ti) 4.0 to 6.0 
and fitting the Ti I line at 6716. 6A. The best x 2 revealed an abundance of log(Ti) = 5.22 
in excellent agreement with the result of Luck & Heiter (2005). The latter value was used 
in the determination of the ^4 ratio. The effects of Ti abundance were investigated further 
by constructing an additional set of models at a given Ti abundance with y|4 ranging from 
0.0 to 0.1. The best fit to the Li line profile provided a ^0 ratio at each Ti abundance. The 
results are shown in Figure 4. The Ti abundance and 7^ ratio are anti-correlated (also see 
Reddy et al. (2002)) in the sense that for a given Li profile, a lower Ti abundance would lead 
to an increased ^4 ratio. However, in the case of GJ 117 studied here, in order to drive the 
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t^t ratio to zero, a Ti abundance of five times greater (log(Ti) = 5.92) would be required. 
We emphasize that both our results and those found in the literature are consistent with 
log(Ti) = 5.22. 

An additional check on any significant blending that may contaminate the red wing of 
the lithium profile was done by comparing the spectrum of GJ 117 to the template star 
HD 131511 which has a similar spectral type. The two objects have similar metallicity 
(Fe/H =0.13 to Fe/H = 0.15) (Luck & Heiter 2005) but HD 131511 has no detectable 
lithium (Takeda & Kawanomoto 2005). HD 131511 was observed at the McDonald 2.7m on 
2005 May 24 with a resolution of ~130,000 and a signal-to-noise just over 120. The spectra 
of the two objects are compared in Fig 5. Features in the 6707.8 to 6708.3 A region were 
examined on the HD 131511 spectrum and showed an equivalent width of less than 2 mA. 
We therefore conclude that our estimate of ^| on GJ 117 is not effected significantly by 
other lines. 

We have also analysed the bisectors of 4 nearby spectral lines of Fe I and Ca I on GJ 117 
and found them to be vertical with an uncertainty of about 0.004 A (0.2 km/s). The red 
asymmetry in the Li profile is therefore not caused by the atmospheric effects of granulation. 
Our findings are consistent with Gray (2005) who analyse the shape of the bisector of the 
K2V star e Eri, a star with atmospheric properties similar to our object. Our findings are 
also consistent with Smith et al. (1998). The lack of asymmetries can be explained by the 
fact that in active stars, like GJ 117, strong magnetic fields inhibit convection which is known 
to be the case in sunspots. 
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2.3. Can the 6 Li abundance be attributed to enhanced activity ? 

Deliyannis & Malaney (1995) have suggested that flare activity may account for 6 Li 
production on HD 84937. However, Lemoine et al. (1997) have questioned the feasibility of 
this scenario on the basis of the energetics involved. We will now examine if the activity 
of GJ 117 can account for the detected 6 Li on this source. As already mentioned in the 
introduction GJ 117 has shown evidence for chromospheric activity. With a rotational ve- 
locity of 5.5 km s _1 GJ 117 is rotating faster than the Sun and has a deeper convection 
zone. Its X-ray luminosity is at least one order of magnitude higher than the X-ray lu- 
minosity of the Sun at solar maximum (Mathioudakis & Doyle 1989; Hempelmann et al. 
1995). Its X-ray emission implies a total average flare energy of f=s 5 x 10 29 ergs s _1 and 
a corresponding total average U-band flare energy of ~ 5 x 10 28 ergs s _1 (Doyle & Butler 
1985; Whitehouse 1985). The most energetic solar flares are the X-class with energies of at 
least 10 32 ergs. The GOES database shows that in 2000, a solar maximum year, 17 X-class 
flares were observed. Assuming that these values also hold for GJ 117 they imply a total 
energy of E« 5 x 10 49 ergs over 1 Gyr. With a production efficiency of 10~ 3 atoms ergs^ 1 
we have 5xl0 46 atoms (Livshits 1997). The mass of the convection zone of GJ 117 is 0.06 
M (Pinsonneault et al. 2001). With an estimated log(Li) = 2.45 and t|4 = 0.03 we obtain 
~ 5xl0 44 atoms. In the largest solar flares the number of protons with energy over 30 MeV 
is more than 10 33 . This value implies a 3xl0 29 7 Li atoms s _1 (Livshits 1997) and 7xl0 28 6 Li 
atoms s^ 1 (Canal et al. 1980) or 2.2 xlO 45 6 Li atoms in 1 Gyr. Our estimates are therefore 
consistent with the 6 Li generation by spallation processes in flare events. 

3. Concluding Remarks 

High resolution VLT UVES observations of the dKl star GJ 117 show the detection of 
6 Li at the 3% level = 0.030±0.010). Our spectral analysis included the Ti I transitions 
and atomic data used by Reddy et al. (2002) for similar studies of the planet hosting star 
HD 82493. GJ 117 is considerably more active than the Sun. We have examined the 
possibility of Li production by spallation in stellar flares and concluded that the energetics 
involved are sufficient for this process to take place. Future work on this subject could include 
the study of the Li profile during flare events where an enhanced 6 Li may be expected. 

"Based on observations collected at the European Southern Observatory, Chile (70. D- 
0075A)". We thank David Lambert for useful discussions and an anonymous referee for sug- 
gested improvements. PHH was supported in part by the Pole Scientifique de Modelisation 
Numerique at ENS-Lyon. Some of the calculations presented here were performed at the 
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Hochstleistungs Rechenzentrum Nord (HLRN), at the National Energy Research Supercom- 
puter Center (NERSC), supported by the U.S. DOE, and at the San Diego Supercomputer 
Center (SDSC), supported by the NSF. DC and MM are also grateful to the Defence Science 
and Technology Laboratory (dstl) for support under the Joint Grants Scheme. 
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Fig. 1.— VLT UVES spectra of GJ 117 showing the Li I 6707.8 A region. 
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Fig. 2. — Ax 2 as a function of the ^ ratio for Li I model line profiles for GJ 117 (see text). 
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Fig. 3. — GJ 117 spectrum for the Li I 6707.8 A region plotted as a solid histogram. Over- 
plotted are PHOENIX models with £g ratios of 0.0 (dotted), 0.03 (dashed), and 0.10 (dash- 
dotted) lines. The lower panel shows x 2 from the difference of the data minus the model 
(see text). 
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Fig. 4. — The Ti abundance as a function of the jj^ ratio (see text). 
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Fig. 5. — Comparison of the VLT GJ 117 spectrum (solid histogram) to that of a K2 template 
star, HD 131511 (dashed line). HD 131511 was observed at the McDonald telescope in May 
2005 (see text). 
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Table 1. Stellar and Lithium Parameters 



Source 


logTe// 


log g 


vsini 

(km s _1 ) 


Li eqw 

(mA) 


log(Li) 


6 Li 


GJ 117 


3.70 


4.6 


5.5±0.2 


211 


2.45 


0.030±0.010 



